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SUMMARY 


An overview of four new analysis capabilities, which have been developed and 
integrated into the GSFC NASTRAN Thermal An^yzer, is given. To broaden the sc<^e 
of applications, these new additions provide, the NTA users with the following capabili- 
ties of: (1) simulating a thermal louver as a means of the passive thermal control, 

(2) simulating a fluid loop for transporting energy as a means of the active thermal 
control, (3) condensing a large-sized finite-element model for an efficient transient 
thermal analysis, and (4) entering multiple boundary condition sets in a single submis- 
sion for execution in steady-state thermal analyses. 


INTRODUCTION 


Four new capabilities have been recently developed and integrated into the GSFC 
NASTRAN Thermal Analyzer (NTA), which is a finite-element based general-pi^pose 
heat transfer computer program in the NASTRAN system (refs. 1, 2, 3). The NTA 
was originally designed to analyze a large and complex space-borne telescope optical 
sub— system (refs. 4, 5). The new additions have broadened the scope of engineering 
application of the NTA from component level to system level, so that complete space- 
craft, which would generally require passive and/or active thermal controls, can be 
treated. Also, the manipulative options relating to the operational efficiency of the 
computer program can become critical when large numbers of elements or grid points 
are employed in modeling rather than dealii^ with a problem of a small size. Specifi- 
cally, the four new capabilities are: 

(1) Temperature-dependent emissivity and absorptivity. 

(2) One-dimensional thermo-fluid elements. 

(3) The condensation of a finite-element thermal model including the non-linear 
radiative effect. 
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(4) The entry of multiple boundary condition sets in a single submission for 
execution. 

The purpose of this paper is to provide NASTRAN users with an overview of these 
new capabilities. Definition, scope, features and limitation of each individual item are 
given. New elements needed to facilitate functional objectives are introduced. The use 
of these new elements together with the appropriate solution routines are described. 
Essential parts of program changes and new elements interfacing with different func- 
tional modules in the NT A are outlined. 


TEMPERATURE-DEPENDENT EMISSIVITY AND ABSORPTIVITY 


In the Level 15- 5 NTA, the solution algorithm pertaining to the radiative exchanges 
is limited to the diffuse-gray case. The programming code is, therefore, sufficient to 
enter a single constant value for the emissivity, e , and the absorptivity, a, of each 
surface. As a result, the radiation matrix (the module RMG) is symmetric. These 
program restrictions have been eliminated by this new capability which permits an 
independent specification of temperature-dependent emissivity e(T) and absorptivity a(T). 
The application of this feature is for transient thermal analysis only. A notable engi- 
neering application is to simulate a thermal louver as a means of the passive thermal 
control. 

To accommodate temperature-dependent emissivity and absorptivity through the 
boundary surface element, HBDY, the RMG module would have to be restructured. 

Such a rigorous approach would result in a prohibitively large effort, and excessively 
long computer run times would be required for the following reasons: 

(1) The matrices which express the radiation properties are generated in the module 
RMG v/hich is executed early in the program. 

(2) If emissivity and absorptivity are temperature-dependent functions, it would be 
necessary to execute the RMG module and all subsequent modules up to and including 
the solution module, TRHT, for each time step. 

(3) If radiative heat flow was requested in the output, it would be necessary to execute 
the data recovery modules after each time step, because NTA does not have the 
ability to store the separate radiation flux matrices [Qge ] for each time step. 

Although less flexible than the rigorous approach, an alternative has been adopted 
to employ the existing general nonlinear capability of NASTRAN (ref. 6). A nonlinear 
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element, not.TNR *. has been newly introduced (ref. 7). It allows the emissivity and 
the absorptivity to be given independently as tabular functions of the radiating surface 
temperature. 

NOLIN5 differs from other NASTRAN nonlinear load cards in that a large number 
of grid points may be involved. An area factor Aj is associated with each grid point. 

For the emitting surface, the convention is adopted that Aj >■ 0. The sum of the area 
factors of the grid points associated with the emitting surface is one; hence the average 
temperature , of the emitting surface is 

T„ • ?A,T, (1) 

Other grid points, assigned to receive the radiation, are given area factors which are 
negative and are not included in the calculation of Tgy . The emissivity e(Tav ) is given 
as a tabular function of Tjv by one of the TABLEDi data cards, as is the absorptivity 
a(Tav ). The rate at which the radiative heat is treated as a nonlinear load at point i is 
given by 


N, = -soA,.cr„)a„ + T.j,)‘> 


( 2 ) 


where S and Ai are data items and a, the Stefan-Boltzmann constant, and T^bs » a datum 
temperature in an absolute scale, are rigid format parameters. The following changes 
have been made to the prc^ram: 

IFP New card NOLIN5 . 

DPD Add NOLIN5 to subroutine DPD4. 

TRHT Add parameter TABS, SIGMA; add NOLIN5 to 

subroutine TRDID. 


Rigid Format (SOL 9) Add parameter to TRHT. 


ONE-DIMENSIONAL THERMO- FLUID ELEMENTS 


A simple but commonly used enei^y transporting means is the fluid loop consisting 
of a flowing fluid inside a tube, which may be made to contact atoermally conductive 
structure. New elements have been developed and added to the NTA to allow the effect 
of an active thermal control to be included in a model. 


♦The names of actual NTA cards are capitalized and underlined. 
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The tube to convey the flowing fluid is modeled using the existing TUBE element 
with modification thnt allows a tube of uniform wall thickness to have a constant bore 
and a linearly tapered geometry as well. The convective heat exchange between the 
fluid and the tube and any heat transferring from the tube to the surroundings are 
modeled using HBDY elements as shown in Figure 1. The HBDY element has been 
modified, as is reflected in the PHBDY data card, to accommodate the linear taper. 

The developed thermo-fluid elements are available to both the cases of steady- 
state and transient thermal analyses. For steady-state analysis, a constant flow rate 
is specified by the user. A new fluid element FTUBE has been introduced. The 
unusual feature of this element is that it possesses an asymmetric conduction matrix 
which is characterized by the fluid phenomenon (ref. 8). This element has created a 
connection C^UBE and a property PFTUBE data cards. Element routines are needed 
for IFP, SMAl, SMA2, SDR2, and OFP. This new fluid element can be used in the 
nonlinear steady-state solution routine, the rigid format APP HEAT, SOL 3. 

For transient analysis, a table, TABLED! , may be provided by the user with a 
time-dependent flow rate function. In addition to the FTUBE element, the fluid cylinder 
is also treated as a nonlinear load. A new nonlinear element N FTUBE has been 
introduced. Both methods can be used in the transient solution routine, the rigid format 
APP HEAT, SOL 9. These two approaches, the element method and the nonlinear load 
method, are equivalent and equally simple to apply as far as the user is concerned. 

The following changes have been made to the program for this thermo-fluid 
capability: 

Modify PTUBE, CHBDY : new CFTUBE. PFTUBE. 
NFTUBE . 

New entries for block data. 

FTUBE routine (as 5 mimetric output); fix HBDY type 
"FTUBE’'. 

Compute convected power in FTUBE 
Output SDR2 results. 

Add NFTUBE to subroutine DPD4 
Add NFTUBE to subroutine TRDID. 


IFP 

GPTAl 
SMAl, SMA2 

SDR2 

OFP 

DPD 

TRHT 
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THE CONDENSATION OF A FINITE- ELEMENT THERMAL MODEL 


Condensation is the reduction of the number of temperature variables (or the 
generalized degrees-of-freedom in the NASTRAN system) used for solution efficiency. 
The OMIT feature of NASTRAN was designed to provide the model condensation. This 
feature has been extended to the transient thermal analysis, using the rigid format APP 
HEAT, SOL 9, t<? include linear and nonlinear boundary conditions. It is accomplished 
through matrix trmsfer and partitioning. The computer code has been modified to 
allow for asymmetric matrix operations. As a result, the restriction used to limit 
"a grid point which is OMITted may not touch any nonlinear elements" has been removed. 

This feature, however, is not available to the nonlinear steady-state thermal 
analysis. Since the possible increase of the computer run time expenditure due to the 
mathematical manipulations would not be justified by the saving that would be realized 
by a condensed model to the steady-state case. 

The program chaises made to accommodate this capability have been modifications 
of MCE2 and SMPl for asymmetric operations. 

THE ENTRY OF MULTIPLE BOUNDARY CONDITION SUBCASES IN ONE EXECUTION 


The counterpart of this capability in NASTRAN has been known to its users as the 
SUBCASE, but this capability was not available to any of the three rigid formats in the 
NTA. To loop through several subcases of disttoct boundary condition sets is a very 
desirable feature, especially for a large-sized NTA model, in which the types of input 
data card of a problem generally consist of the conduction elements (the connection and 
property data cards), the boundary surface elements, the material cards, the radiation 
matrix input cards (RADLST and RADMTX ), and different thermal loading cards for 
specifying boundary conditions. Usually, thermal loadings need only a few cards to 
describe needed quantities and conditions in steady-state cases, while the rest of cards 
constitutes the major part of a bulky input deck. This feature has been developed and 
used to eliminate time-consuming repetitive submissions of bulky input data decks in 
succession with only changes pertaining to thermal loadings varied from case to case. 
This multiple entry feature is limited, however, to steady-state thermal analyses using 
the rigid formats of APP HEAT, SOL 1 and SOL 3. Specifically, a DMAP- ALTER has 
been added to handle the following loops for 

(1) Only thermal loads are changed. 

(2) The MPC and SPC constraints are changed. 

(3) The boundary surface elements CHBDY are changed. 
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The required program changes for looping in this rigid format APP HEAT, SOL 3 are 
summarized as follows: 

Write ALTER 

Fix to IFPl, PARAML for ALTER. 

The DMAP- ALTER packet is shown in Figure 2 for reference. 
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Figure 1. Combination of the NT A elements to simulate thermo-fluid effects 
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s THE DMAP-ALTER PACKET TO ALLOW BOUNDARY CONDITION S^TS CHANGED IN RF-hEAT ,3 

s TO ALLOW THERMAL LOAD TO VARY FROM SUBCASF TO ^SUHCrSb^^ALTERSa 
ALTER 7,7 

EOUIV j?|^^^PECT^/NHBOY/HSIL,PSlL/NHBOY/HEOEXIN,PEOIN/NHBDY/HBGPnT, 
ALTER13,13 

t!^i-IPA'^»HGPSETS,HELSETS,CASECC*PBr,Pf)r,PFoi.j.psil / 

ALTER 19 ^'^'-^^^^^^^^^'-/^^'-USET/JUMPPLOT/V,N,HPLTFLG/V,NWHPFiLE^f 

SAVE AAPEATTfS§i:&*i’'*”‘''''''''’"’“'=''EAIT/V,N,.vnLW)P S 
PARAML CASEXX//DT I/l/ 8 //V,N,TFMPMATF f. 

ALTER 3^^|^S«/V’^'*TEMPMATE;-10 * 

GP4 CA^E XX ,GE 0M4 , HEOtX IN ,HS I L »H 6 PDT . • /HRG. .MU';-? /hihcut/ 

ALTEA 

ALtIr^3?‘;37’'^'^‘'''''^’'^^^'^^^^’^0'^J^’^'^‘:''CT,HN0SET,M:.|()L,hN0A S 

ALTER 59 » ^^^0/ V , Y , KSYM/- 1 /O $ 

^^SLT,H8GP0T,HCSTM,HSIL,HEST ,MPT,HGPiT,i:Dr, , lAScXX, DI T/hpg/ 

ALTER 63,63 

ALTER ’’HI'G/ttHPS.HPF $ 

SASEXX,HCSTM,MPT,DIT,HE0EXIN,,HGPTT,FDT.HHGPI)P HPr Hi,r 
ALTER 0 ^ ^ ^ ^ 1 HOOGl , HOUGV 1 , HOr S I , HO 1 , J G V I / S f ,A T fc S % 

RLTSET PC06 r HEOEX IN ,HEC T/ PSMtS ♦ DP I Th> A(^ . nnw> «: T c . i^v^i c L. T c y , 

V , N , ( J J $ ■ ' ~ - 7 i_ , y V , -V , 1 / 

SAVE DSIL $ 

SAVE HPhlLri^''^^'^^ * 

ALTER 83 

COND fin IS, REPEAT! 

REPT L00PT0P,100 
JUMP HERR0R3 $ 

$ END OF DMAP-ALTER (53$48) 


Figure 2. The DMAP-ALTER packet for the multiple boundary condition subcases ir 
one execution. 
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